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Jh ■ Abstract. We study the Josephson oscillations of two coupled elongated condensates. Linearized calcula- 

tions show that the oscillating mode uniform over the length of the condensates (uniform Josephson mode) 
is unstable : modes of non zero longitudinal momentum grow exponentially. In the limit of strong atom 
interactions, we give scaling laws for the instability time constant and unstable wave vectors. Beyond the 
^ ' linearized approach, numerical calculations show a damped recurrence behavior : the energy in the Joseph- 

, son mode presents damped oscillations. Finally, we derive conditions on the confinement of the condensates 

• . to prevent instabilities. 
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O • vortices and topological excitations - 03.75.Kk Dynamic properties of condensates; collective and hydro- 
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^ ; 1 Introduction 

*n : 

i Josephson oscillations arise between two Bosc-Einstcin condensates coupled by tunneling effect. They have been ob- 
£N| i served in superfluid HeliumpQ and in superconductors and have recently been achieved in dilute atomic BEC in a 
<0 double well potential^- The physics of two coupled condensates has been extensively studied in a two modes model, 
in where only two single particle modes are involved For atoms interacting in each well through a two-body inter- 
■ action, different regimes are reached depending on the ratio between the tunneling strength to the interaction energy 
c/3 | of atoms in each well[HE]. For small interaction energy, one expects to observe Rabi oscillations. For large interaction 
. — ' energy one enters the Josephson regime. In this regime, oscillations around equilibrium configuration have a reduced 
c/2 ' amplitude in atom number and their frequency depends on the mean field energy. Finally, for very large interaction 
energy, quantum fluctuations are no longer negligible : the system is in the so called Fock regime and oscillations 
of atoms between the wells do not occur any more. In this paper, we assume this regime is not reached. Oscillations 
between the two wells, both in the Rabi and in the Josephson regime, are then well described by a mean field approach. 

Atom chips [7j are probably good candidates to realize Josephson oscillations of Bose-Einstein Condensates as they 
enable the realization of micro-traps with strong confinement and flexible geometries. A possible configuration to 
realize a tunnel coupling between BEC on an atom-chip is proposed in In this proposal, the two condensates are 
1 very elongated and are coupled all along their longitudinal extension. With such an elongated geometry, both the 
Rabi and the Josephson regime could be accessed. However, in this case, tunnel coupling may be larger than the 
longitudinal frequency and the two modes model a priori breaks down. In this paper, we are interested in the stability 
of the uniform Josephson mode where all the atoms oscillate between the two wells independently of their longitudinal 
position. In the absence of interaction between atoms and if the transverse and longitudinal trapping potentials are 
separable, the longitudinal and transverse degree of freedom are decoupled and one expects to observe stable Rabi 
oscillations between the condensates. On the other hand interactions between atoms introduce non linearities that may 
couple the two motions. For a homogeneous situation such as atoms trapped in a box-like potential, uniform Josephson 
oscillations are a solution of the mean field evolution equations and are a priori possible, even in presence of interactions 
between atoms. However, the non linearities introduced by interactions between atoms may cause instability of this 
uniform Josephson mode. Similar modulational instabilities appear in many situations of nonlinear physics such as 
water waves propagation^] or light propagation in a non linear fiber|llj. In the context of Bose Einstein condensates, 
they have been observed in presence of a periodic potential, at positions in the Brillouin zone where the effective mass 
is negative |12U13lll4j . In our case a modulational instability would cause uniform Josephson oscillations to decay into 
modes of non vanishing longitudinal momentum. The goal of this paper is to investigate those instabilities. 
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We assume that all the relevant frequencies (interaction energy and tunnel coupling) arc much smaller than the 
transverse oscillation frequencies in each well so that we can consider only a one dimensional problem. Thus, the 
system we consider is described by the Hamiltonian 



H 



+U(z) [t/>l(z)tl>i{z)+il>l(z)M*) 
ipl(z)Mz) + 4(z)Mz)]}, 



"7 



(1) 



where g is the one-dimensional coupling constant and U(z) is the longitudinal trapping potential. For a harmonic 
transverse confinement for which lo± <C fr 2 /(ma 2 ), we have g = 2fiw_^a, where a is the scattering length [T5]. The 
parameter 7 describes the tunnel coupling. 

We are interested in the stability of uniform Josephson oscillations around the equilibrium configuration where the 
two condensates have the same phase and equal longitudinal density. In the sections 2-4, we consider a homogeneous 
configuration where U(z) = 0. In the sections 2 and 3, we calculate the linearized evolution of modes of non zero 
longitudinal momentum in the presence of uniform Josephson oscillations. In the section 2, we give results of a 
calculation valid both in the Josephson and in the Rabi regime. In section 3, we show that, in the Josephson regime, 
the system is well described by a modified Sine-Gordon equation. For small amplitude oscillations, we derive scaling 
laws for the instability time constant and the wave vectors of the growing modes. In section 4, we go beyond the previous 
linearized approaches and present numerical results. We observe damped oscillations of the uniform Josephson mode 
amplitude. Such oscillations are similar to the Fermi-Pasta-Ulam recurrence behavior |lf>lll7j . In the last section (5), we 
present numerical calculations in the case of a harmonic longitudinal confinement. We show that Josephson oscillations 
are stable for a sufficiently strong confinement and we give an approximate condition of stability. 



2 Numerical linearized calculation 

To investigate whether Josephson oscillations arc unstable with respect to longitudinal excitations, we use a linearized 
calculation around the time-dependent solution corresponding to uniform Josephson oscillations. Writing ipi = ipi+S^i 
and "02 = ^2+^2 with uniform (z-independent) ^1,2, Eq^gives to zeroth order the coupled Gross-Pitaevski equations 

ift^rfi = 9\<Pi\ 2 <Pi - lf2 + (7 - PosVi , 

f • ( 2 ) 

ih 'dl, (P2 = 9 ^ 2(p2 ~ 1Vl + (7 ~ P05V2 

We shifted the zero of energy by adding to the Hamiltonian a "chemical potential" term "f — pog, where po is the density 
of each condensate at equilibrium. We recover here the well known results established for a two modes model|4*||^]|(J). 
More precisely, writing (f \ = N\jh e 101 and <p% = \J (N — N\)/L e 102 where L is the size of the system, Eq|21 implies 
that the conjugate variables d\ — 62 and k = (Ni — A^)/2 evolve according to the non rigid pendulum Hamiltonian 
Hp = Eck 2 /2 + Ejy/l — Ak 2 /N 2 cos(#i — #2) where the charge energy is E c = 2g/L and the Josephson energy is 
Ej = 7N. We consider oscillations of 9\ — 62 around of amplitude OSC - Let us now consider the evolution of 
excitations around those uniform oscillations. To first order in 0^1,2, EqQ]yiclds the coupled Bogoliubov equations 



(3) 





where, for i = 1,2, 



2 



and the coupling term is 



d=\ -ot + 2 ^I 2 -pos + 7 -gn (4) 

9<p? hih- 2 9\^\ 2 + p»g-i) 

<-(7°)- m 

Instabilities arise if there exist modes growing exponentially in time under Eq|3J The evolution matrix is invariant 
under translation so that we can study independently plane waves modes e lkz (ui, Vi, U2, V2), the second derivatives 
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Fig. 1. Evolution of the square amplitude of the symmetric (fat lines) and antisymmetric (thin lines) excitations of wave vector 
k — 0.1^/mpog/h (a) and k — 0.3^/mpog/h (b). Those graphs are computed for 7 = O.lpoS and a uniform Josephson oscillation 
amplitude OS c = 0.6. 



in Ci and C 2 being replaced by — k 2 . Note that the evolution of excitations depends only on the four parameters k, 
pog, 7 and OSC . For a given k component, we numerically evolve equations |2] and |3 FigQ gives the evolution of the 
square amplitude of the symmetric mode |m s | 2 = |tii + U2I 2 and of the antisymmetric mode |u a | 2 = |ui — W2I 2 for 
two different k vectors, for 7 = 0.1po<? and for <9 osc = 0.6. For these calculations, we choose the initial condition as 
(ui,vi, 1*2, K2) = (1, — 1) — 1) 1)- In the two cases, wc observe a fast oscillation at a frequency close to the frequency 
of the antisymmetric mode y/ {2pog + 27 + h 2 k 2 /2m)(2~f + H 2 k 2 /2m) and a slower oscillation at a frequency close to 
that of the symmetric mode -v / (2pna + h 2 k 2 I2m\h 2 k 2 /2m|18|. On top of this, we observe, for k = 0.1, an exponential 
growth e 2Ft of |ui +M2I 2 and \u± — M2I 2 , signature of an unstability. We find that, for given pog and (9 osc , the instability 
domain in k is [0, k max \. FigEl gives the maximum growth rate r and the maximum unstable wave vector k max . 



3 Calculation in the Josephson limit 

In this section, we focus on the Josephson regime where 7 <IC ono|10|. In this regime the amplitude of oscillations in the 
relative density Sp remains small compared to the mean density and one can assume p\ = P2 in the Josephson energy 
term of the Hamiltonian. Furthermore, wc restrict ourselves to long wavelength excitations described by phonons and 
we neglect anharmonicity of phonons. Then, the Hamiltonian reduces to 

Hj = H s + H SG + H c , (6) 

where, writing i^i = ip 2 = v^e^ 2 , Q a = Q\ - 6> 2 , S = 0\ + 9 2 , p a = (pi - Pi) 12 and p s + p = (p 1 + p 2 )/2, 

describes the symmetric phonons, 
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Fig. 2. Maximum instability rate of excitations (a) and maximum wave vector k of unstable modes (b) as a function of the 
amplitude of the relative phase oscillations for 7 = 0.05poff (stars and solid line) 7 = O.lpofl (crosses and dashed line) and 
7 = 0.2po3 (circles and dotted line). The points are the results of the linearized numerical calculations presented in section^ 
and are given with a precision of 10%. The continuous lines are given by diagonalising the four by four matrix as presented in 
section^ 



is the Sine-Gordon Hamiltonian and 

H c = -2 7 J Ps {cos(6 a ) - l)dz (9) 

is a coupling between the symmetric and antisymmetric modes. The Sine-Gordon Hamiltonian has already been intro- 
duced in the physics of elongated supraconducting Josephson junction|2j. In those systems, symmetric modes would 
have a very large charge and magnetic energy and do not contribute. The Sine-Gordon model has been extensively 
studied [T5|. In particular, it has been shown that, for a Sine-Gordon Hamiltonian, oscillations of well defined momen- 
tum (in particular k = 0) present Benjamin- Fcir instabilities |19). Our system is not described by the Sine-Gordon 
Hamiltonian because of the presence of H c . In the following, we derive results about stability of our modified Sine- 
Gordon system. As we will see later, we recover results close to that obtained for the Sine-Gordon model. 

The Josephson oscillations correspond to oscillations where p a = p osc and 8 a = 9 osc are independent of z. They 
are given by 

27P0 sin(0 

osc) tl 

(10) 



dt 

de ot 



at 



They also induce an oscillation 60% of 6 S given by 



86, 



is) 



dt 



-2 7 (cos(0 osc )-l)/fi. (11) 



To investigate whether some non vanishing k modes are unstable in presence of a Josephson oscillation, we linearize, as 
in the previous section, the equation of motion derived from EqEJaround the solution p OSCl OSC . Because of translational 
invariance, we can study independently the evolution of modes of well defined longitudinal wave vector k. Writing 
Pi = Po + Posc + (Spa + 5p s )e ik *, P 2 = Po- Pose + (Sp a + Sp s )e ikz , 9 1 = {9ot + 6 osc + (59 s + S6 a )e ikz )/2, and 
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J osc 



S9 a )e )/2, we find the evolution equation 



dt 



( Sp a /po\ 

se a 

Sps/pO 

\se s 

I + 2tcos(6» os 

-2poS 


27sin(6l osc ) 



7 



27 sin(Sosc) 



-2poff 



-',.2 



2m 





/ & pal pa \ 
S6 a 

Sps/po 

\se a ) 



(12) 



We solved numerically Eq^J and Eq^| and we find that modes of low k wave vectors are unstable. Fig|3] gives the 
instability rate and the maximum k wave vector of unstable modes. Those results agree within 10% to the more general 
results of the previous section as long as < 0.2 and the oscillation amplitude fulfills osc < 0.6. 

To get more insight into the physics involved and to obtain scaling laws for the instability rate and the instability 
range in k 7 we will perform several approximations. The evolution matrix M of Ea ll2l is periodic in time with a period 
uj. We can thus use a Floquet analysis PH] and look for solutions of Eq^]in the form 



+00 

E 



+00 



& inwjt c _ e ivt 



~y ' ginuijt 



n— — 00 



n— — 00 




Expanding Eg 1 121 for each Fourier component, we find 



= —Lo,jnc n — iMoCn — i M. 



m^n—m: 



where the time independent matrices M n are the Fourier components 



M„ 



2nh 



-imujt 



M(t)dt 



(13) 



(14) 



(15) 



Thus, solutions of Eq^|are found as eigenvalues of the linear set of equations (|14f) . The mode is unstable if there 
exists an eigenvalue of non vanishing real part and its growth rate is the real part of the eigenvalue. 

For OSC = 0, only the dc component M is not vanishing and its eigenvalues are ±iv a = ±iy / 2p g(2^ + hk 2 /2m) 
and zLlu s = ztihk^ pag/m corresponding, for each Fourier component n, to the symmetric modes and antisym- 
metric modes c± ■ The four states c_ _ 1 , c_ Q , c_|f ^ and c+ x form a subspace almost degenerate in energy and of 
energy far away from the other states as depicted Fig@J Thus, we will restrict ourselves to those states in the following. 
In the limit of oscillations of small amplitude OSC , the matrix elements of M can be expanded to second order in 
9 osc . Furthermore, the oscillations are well described by OSC = OSC cos(wji), where ujj = 2y/jp g(l — 0l sc /16)/h. We 

then find that, in the 4 dimensional subspace spanned by (c_ _ 1 , c_ ' Q , c^!f ' c +\)' ^ nc c ig enva hte v of Eq^]arc the 
eigenvalues of the four by four matrix 

I i(-u a + wj + ^01^/4) -i 1 osc f a /{2f s ) 1 osc f a /{2f s ) -J0 2 osc f 2 a /& \ 

hOoscfa/fs/2 IUJ S -J0 OS cfa/(2f s ) 

-f0oscfa/(2fs) -i<jj s -i'J®oscfa,/('2fs) 

-lOl sc fl/% 10oscfa/(2f s ) tJ0oscfa/(2fs) »(w„ - Uj - 7^c/a/4) J 



M 



(16) 



where f a = (2p g/ (h 2 k 2 /2m + 2 7 ))( x / 4 ) and f s = {4mp g/h 2 k 2 )^/ 4 \ 

We numerically diagonalise this matrix and find the instability rate as the largest real part of the eigenvalues. For 
a given oscillation amplitude OSC , scanning the wave vector k, we find the largest instability rate and the maximum 
wave vector of unstable modes. FigOU compare those results with the values obtained by integration of EqEl We find 
a very good agreement in the range 9 < 0.6 and 7/pog < 0.1. Finally, in FigEl we compare the instability rate and the 
maximum unstable wave vector found with this simplified Floquet analysis with the more general results of section 1. 
We find a very good agreement as long as 7/poff < 0.2 and <9 osc < 0.6. 

If we restrict ourselves to terms linear in <9 osc , then the only effect of the Josephson oscillations is to introduce a 
coupling between the symmetric and antisymmetric mode. We checked that this coupling alone does not introduce any 
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Fig. 3. Comparison between numerical evolution of Eas YTffl and XT^ (points) and the results obtained by diagonalising the 4 by 
4 matrix of the Floquet representation (lines). Parameters are 7 = 0.1 x pog (stars and continuous lines) and 7 = 0.05 x pog 
(crosses and dashed line). 




n = -1 n — n = 1 

Fig. 4. Floquet representation of the equation Ea MZA The ellipse surrounds the four states that are considered in the calculation 
of instability rates. 



instability. Thus instability is due to the quadratic terms. Those terms contain a modulation at 2u)j. This modulation 
corresponds to the modulation of the frequency of the antisymmetric mode 

u£ = 2 Po g(h 2 k 2 /2m + 2 7 - 2-/0 2 osc /4) + 7Po<7©L cos(2lj jt). (17) 

This parametric oscillation leads to instability for k <G [0, osc y / mj/2/h] and the instability time constant at resonance 

is r = 2 sc ^-^^- ■ We recover here the well known results of Benjamin- Feir instability derived for example in JHj using 
the multiple-scale perturbation technique. In our case, the coupling to the symmetric mode will modify those values. 
However, for small values of 7, the qualitative behavior is unchanged. Indeed, as seen in FigEl as long as 7 < 0.05po.9 
and within a precision of 10%, the instability rate _T scales as 

r = O.122(1)0 O 2 SC V™A (18) 

and the maximum wave vector of unstable modes as 

k max = O.97(l)^p0 osc . (19) 
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Fig. 5. Maximum instability rate normalized to the Josephson oscillation frequency (a) and maximum wave vector of unstable 
modes normalized to ^/mrj/h (b) as a function of the oscillations amplitude Oosc for different ratios "f/pog (from lower curves 
to upper curves : 0.02,0.06,0.1,0.14). Fat dashed lines are the scaling laws Eq |18l and 1191 Thin continuous lines are found by 
diagonalising the matrix of Ea llfil 

For larger 7, the r and k max are higher than those lows as seen in FigEl 



4 Beyond the linearisation 

The two previous sections give a linearized analysis of the evolution of perturbations. They show that the presence 
of uniform Josephson oscillations produces instabilities of modes of non vanishing momentum. The energy in these 
mode grows and consequently, the energy of the uniform Josephson mode decreases and one expects a decrease of the 
uniform Josephson oscillations amplitude. Such a decrease is beyond the previous linearized analysis and we perform 
full numerical calculation of the evolution of the mean fields ip\{z,t) and ^2(2, t). The evolution equations derived 
from Eq^are 

d , h 2 d 2 ^ 



ih dt^ 



2m dz 2 
h 2 d 2 ip 2 



fflV'ilVi 



7^2 



(20) 



2m dz 2 

Fig E| gives the evolution of the total number of atoms in the condensate 1, Ni = j\ipi\ 2 , for initial amplitude 
osc = 0.6 and for different values of 7/(poff)- F° r these calculations, the initial state consists in a z-independent phase 
difference OSC between ipi and tp2 superposed on thermal fluctuations of the density and phase of the two condensates 
corresponding to a temperature kgT = p g /10. Wc observe that the amplitude of the Uniform Josephson Oscillations 
presents damped oscillations. For 7 <C p g, the period of these amplitude oscillations is about three times the inverse 
of the instability rate of Eq^| The ratio between the Josephson frequency and the frequency of these amplitude 
oscillations is about 20 and is almost independent on the ratio between 7 and p$g as long as 7 < p§g. For larger 7, 
this ratio increases and more Josephson oscillations are seen in a period of the amplitude modulation. Such amplitude 
oscillations are a reminiscence of the Fermi-Ulam-Pasta recurrence behavior observed in many non linear systems with 
modulational instabilities j 1 911 1 611 1 7) . In particular, this recurrence behavior has been seen in numerical evolution of the 
Sine-Gordon Hamiltonian|21| . In our case, we observe an additional damping which results probably from the coupling 
to symmetric modes. 

The case of an initial amplitude <9 osc = 7r/2 is of particular interest as, in absence of interactions between atoms, 
it corresponds to Rabi oscillations of maximum amplitude. Fig0 gives the evolution of N\ for 7 = pog and an initial 
amplitude OSC = tt/2. 
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Fig. 6. Evolution of the number of atoms in the condensate 1, normalized to the total number of atoms, as a function of 
time. The initial state corresponds to a phase difference between the condensate OS c = 0.6 superimposed on phase and density 
fluctuations corresponding to a thermal equilibrium at temperature ksT = O.lpog. For this calculation, 7 = O.lpoS (a) o,nd 
7 = Pog (b). 
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Fig. 7. Evolution of the number of atoms in the condensate 1, normalized to the total number of atoms, as a function of 
time for 7 = po<? and an initial phase difference between condensates OSC = ir/2. Initial thermal population of excited modes 
corresponding to ksT = O.lpoS is assumed. 



5 Case of a confined system 

In the previous sections, we considered large and homogeneous systems. We found that unstable excited modes are 
those of low wave vectors. In the Josephson limit where 7 -C po9, we derived the scaling law Ea 1 191 for the maximum 
unstable wave vector. In a cloud trapped in a box like potential of extension L, the minimum k value of the excitation 
modes is 2tt/L. Thus, if 

2nH , . 

L < , (21) 

1.0y/rrvy& osc 

the minimum wave vector of excited modes is larger than the maximum unstable k value Ea ll9l and the system is stable. 
This condition can be understood in a different way : the energy of the lowest longitudinal mode is h2TTy/p g/(mL) ( 
here we assume L> hj y/mpog). Thus, we find that the system is stable provided that the energy of the lowest excited 
mode satisfies E exc > O.52luj0 osc where ujj = 2y/jpog/h is the Josephson frequency. 
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Fig. 8. Josephson oscillations of clouds trapped in a harmonic potential of frequency ui — O.lpog/h (solid line) and to = pog/h 
(dashed line), where po is the peak linear density in each condensate. The initial phase difference between the condensates is 
7r/2 and the tunnel coupling is 7 = pog. Ni is the number of atoms in the condensate 1 and N to t the total number of atoms. 

An approximate condition of stability of Josephson oscillations in the case of a cloud trapped in a harmonic 
longitudinal potential of frequency u> is found as follows. The size of cloud, described by a Thomas Fermi profile, is 
L = 2fj,/(moj 2 ), where \i = pog is the chemical potential and po the peak linear density. Then, from the same argument 
as above, one expects to observe stable oscillations for 

lo > cty/jp g© osc = a0 osc toj/2 (22) 

where a is a numerical factor close to one. We performed numerical simulations of the evolution in the case of a 
harmonic potential, adding to both left hand sides of EqsEOI a trapping potential l/2muj 2 z 2 . The initial situation 
is the Thomas Fermi profile superposed on thermal random fluctuations and a global phase difference between the 
condensates O osc = ir/2. The tunnel coupling is 7 = pog. The resulting Josephson oscillations are shown in FigEI 
for co — pog/h and 10 = Q.lpog/h. We observe that for lo = pog/h, Josephson oscillations are stable whereas, for 
lo = 0.1pog/h, oscillations are unstable. 

6 Conclusion and prospects 

We have shown that Josephson oscillations of two coupled elongated condensates are unstable with respect to ex- 
citations of longitudinal modes. The unstable modes are those of small wave vectors. In the Josephson limit where 

7 <C pog, we have derived the scaling lows Eq^Jand EqEHfor the instability time constant and wave vectors. Since the 
frequency of Josephson oscillations are 2^fph~g 1 the first equation tells us that the number of oscillations that can be 
observed scales as 2 sc and is independent on ^y/p^g. This is true as long as 7 < pog. For larger j/^pog), the Josephson 
condition is not fulfilled. Effect of interactions is less pronounced and more oscillations can be observed. Performing 
numerical calculations beyond the linearized approach, we have shown that the system presents a recurrence behavior, 
although it is damped quickly. Finally, we investigated the stability of oscillations in finite size systems. Eg 1211 gives 
the maximum longitudinal size of confined condensate that enables the presence of stable Josephson oscillations. We 
also considered the case of harmonically trapped cloud and give an approximate condition on the oscillation frequency 
to have stable Josephson oscillations. 

The results of this paper are not changed drastically for finite temperature. Indeed, although elongated Bose- 
Einstcin condensates present thermally excited longitudinal phase fluctuations [22E3 , it has been shown in ^2] that, 
because the antisymmetric modes present an energy gap, thermal fluctuations of the relative phase between elongated 
coupled condensates are strongly suppressed. 

Among the possible extensions of this work, two questions are of immediate experimental interest. First, the 
effect of a random longitudinal potential could be investigated. Indeed, it has been proposed to realized elongated 
coupled condensates using magnetic trapped formed by micro-fabricated wires[H], but, for such systems, a roughness 
of the longitudinal potential has been observed 24 25,20 . If the amplitude of the roughness potential is smaller than 
the chemical potential of the condensate, one expects to still have a two single elongated condensate. However, the 
roughness of the potential may change significantly the results of this paper. Second, the effect of correlations between 
atoms may be studied. Indeed, for large interactions between atoms, correlations between atoms become important. 
More precisely, for po < mg/H 2 , a mean field approach is wrong and the gas is close to the Tonks-Girardeau regime |2"?1 
I28U29| . Such a situation is not described in this paper in which a mean field approach has been assumed. Thus, a new 
study should be devoted to the physics of coupled elongated Tonks gas. 

Dynamical instabilities of the uniform Josephson mode are not the only effect of non linearities in the system of 
two coupled elongated condensates and other interesting phenomena are expected. For instance, ref.[3U| shows that 
Josephson vortices similar to the solitons of the Sine-Gordon model exist for large enough interaction energy. 
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